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Abstract We have developed a radiative transfer code, 
CMFGEN, which allows us to model the spectra of mas- 
sive stars and supernovae. Using CMFGEN we can derive 
fundamental parameters such as effective temperatures 
and surface gravities, derive abundances, and place con- 
straints on stellar wind properties. The last of these is 
important since all massive stars are losing mass via a 
stellar wind that is driven from the star by radiation 
pressure, and this mass loss can substantially influence 
the spectral appearance and evolution of the star. Re- 
cently we have extended cmfgen to allow us to un- 
dertake time-dependent radiative transfer calculations 
of supernovae. Such calculations will be used to place 
constraints on the supernova progenitor, to place con- 
straints on the supernova explosion and nucleosynthe- 
sis, and to derive distances using a physical approach 
called the "Expanding Photosphere Method". We de- 
scribe the assumptions underlying the code and the 
atomic processes involved. A crucial ingredient in the 
code is the atomic data. For the modeling we require ac- 
curate transition wavelengths, oscillator strengths, pho- 
toionization cross-sections, collision strengths, autoion- 
ization rates, and charge exchange rates for virtually 
all species up to, and including, cobalt. Presently, the 
available atomic data varies substantially in both quan- 
tity and quality. 
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1 Introduction 

Massive stars are a crucial ingredient of galaxies, and 
the universe. They enrich the interstellar medium 
(ISM) with metals, either through "quasi-steady" mass 
loss, or when they explode as supernovae (SNe). They 
deposit momentum and energy into the ISM, and ionize 
the surrounding gas, producing colorful nebulae (Hn 
regions) . Massive stars also exhibit a wide range of in- 
teresting phenomena including radiation driven winds, 
and colliding winds. Colliding winds, which occur in 
massive star binaries, generate strong shocks t hat give 
rise t o hard X-ray emission (see, e.g., review by Stevens! 
2005). The key tool for understanding massive stars is 
spectroscopic analysis. 

Using spectroscopic analysis we generally wish to in- 
fer the fundamental parameters that describe the star 
— its mass, effective temperature, radius, surface grav- 
ity, and surface abundances. From these, and with the 
aid of evolutionary models, we try to infer both the 
previous evolutionary history of the star and its future 
evolution. As we generally can't study the evolution of 
a single star we need to infer/constrain the effects of 
complex physical processes by studying groups of stars. 
Constraints provided by these studies can then be used 
to improve evolutionary calculations. 

There are still significant uncertainties in evolution- 
ary models (e.g., treatment of convection, stellar rota- 
tion, and time-dependent mass loss). Rotation, for ex- 
ample, is now recognized as a crucial factor in massive 
star evolution — it affects the evolutionary lifetime of 
the star, surface abundances, and the star's evolution, 
and there is a strong interaction between rotation and 
mass loss (e.g., Mevnet and Maeder 2000h . Since stars 
form with a range of rotation rates, it is no longer pos- 
sible to ass ign a position in the HR-dia gram to a unique 
mass (e.g., Mevnet and Maeder 2000f ). 

There are also uncertainties in spectroscopic analy- 
ses. Some of these arise from model assumptions and 
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inadequacies in the atomic data. Others arise from an 
incomplete understanding of the stars we are trying to 
model. O stars, for example, show strong evidence for 
microturbulent velocities approaching the sound speed, 
while Of stars show evidence for macroturbulent ve- 
locities in excess of the sound speed. The origin of 
these velocity fields is unknown. Mass-loss rates are a 
key ingredient of stellar evolution models, but deriving 
accurate mass-loss rates is difficult. While wind the- 
ory provides a qualitative understanding of mass loss 
in O stars, there are fundamental uncertainties since 
theory predicts, and observations show, that the winds 
are highly inhomogeneous. As a consequence of the in- 
homogeneities, mass-loss rates derived from observation 
depend on the diagnostic used. For more information 
on radiation driven winds, mass-loss rates, and prob- 
le ms modeling O stars, the reader is referred to rev iews 
bvlHillierl liooi h lPuls. Vink and Naiarrol |2008f ) and 
Owockil foOoA 



2 CFMGEN — A Spectroscopic Tool 

In order to facilitate analysis of hot stars (stars in which 
molecules and energy transport by convection at the 
stellar surface can be neglected) we have develo ped a 
non-L TE radiative tra nsfer code, CMFGEN (Hilfier ll987l ; 
1990l : Hillier & Miller |l998|) . 

The primary purposes of CMFGEN are as follows: 

1. To derive accurate stellar parameters and abun- 
dances for comparison with "evolution" calculations. 

2. To provide accurate EUV (i.e., A < 912 A) radiation 
fields for input to nebular photoionization calcula- 
tions. 

3. To provide fundamental data for the study of star- 
bursts, star formation in galaxies, etc. 

4. To provide a better understanding of the hydrody- 
namics of stellar winds. 

5. To provide distances to Type II SNe using the ex- 
panding photosphere meth od (EPM — e.g.. East- 
man, Schmidt, and Kirshner 19961 ; Baron et «Lll2004 ; 



Dessart and Hillierl 1 2006T ) and its variants 



6. To provide diagnostics of SNe which can place con- 
straints on the progenitor and the explosion. 

7. To allow the development and testing of approxi- 
mate methods that can be used in more complex 
geometries and in inhomogeneous media. 

cmfgen has been used to study O Stars; Wolf-Rayet 
(W-R) Stars; Luminous Blue Variables (LBVs); A & B 
supergiants; and Type I and Type II Supernovae. 

cmfgen was originally designed to model spectra 
of W-R stars which have a dense stellar wind. The 



hot W-R core, together with the dense stellar wind, 
gives rise to an optical spectrum dominated by emis- 
sion lines — the antipathy of normal spectra in which 
the optical spectral region is dominated by absorption 
lines. Because of the stellar winds, with flow velocities 
of order 1000 km s -1 , it is convenient to solve the radia- 
tive transfer equation in the comoving frame — in this 
frame the opacities and emissivities can be assumed to 
be isotropic. 

In cmfgen, the primary radiation transport equa- 
tions to be solved (assuming V oc r and spherical ge- 
ometry) are 
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are moments of the radiation field, D/Dt is the La- 
grangian derivative, I v is the specific intensity at time 
t, location r, frequency v, and in direction fj,, fj, = cos 9 
where 9 is the angle between the radius vector at r 
and the specific intensity, \ v is the (frequenc y depen- 
dent) opacity, and r\ v is the emis sivity (e.g., iMihalas 
19781 : iMihalas and Mihalasl Il984h . When written in 
this form, the equations are deceptively simple — in 
reality they are very complex. First, these equations 
need to be solved at a large number of frequencies, typ- 
ically 100,000. Due to the djdv term, these equations 
are explicitly coupled in frequency. Second, much of 
the physics is hidden in \v an d f\v In the best case 
scenario these are determined by the local temperature 
and density, but even then the temperature is deter- 
mined by the radiation transport, and thus there is an 
implicit coupling of the opacities and emissivities with 
the radiation field. 

In stellar atmosphere modeling we can neglect the 
Lagrangian derivative term, and the equations are eas- 
ily solved (for non-Hubble flows additional terms are 
introduced). However, for Type II SNe, it has be- 
come increasingly apparent that to accurately model 
the spectra the full time-dependent radiation transport 
equations must be solved. 
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3 LTE & non-LTE 

A great simplification, applicable to most main-sequence 
spectral types, is the assumption of local thermody- 
namic equilibrium (LTE). With this assumption we as- 
sume that the state of the gas is entirely determined by 
its temperature and density (and composition). Thus 
we can use statistical arguments to determine the ion- 
ization state of the gas (Saha equation) and level popu- 
lations (Boltzmann or Saha-Boltzmann equation). The 
LTE assumption is valid when collisional processes 
(which couple directly with the local gas) dominate 
over radiative processes (which couple the gas to gas 
elsewhere by radiation). 

In LTE, we require line lists and opacities (due to 
bound- free, free- free, and bound-bound processes). The 
opacities, for a given composition, are simply functions 
of the density and electron temperature. Deep inside 
the star radiation diffuses, and the transport of radia- 
tion is dictated by the Rosseland mean opacity defined 

by 
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In the photosphere the situation is very different — 
use of the Rosseland mean opacity is no longer appli- 
cable and we must solve the radiative transfer equation 
at every frequency. To do this we need a detailed de- 
scription of the opacities (Fig. [TJ. Since the theoreti- 
cal spectrum is to be compared with observation, it is 
important the bound-bound transitions have accurate 
wavelengths. 

In O stars, W-R stars and SNe (as well as many 
other objects) the assumption of LTE is invalid. In- 
stead we are forced to solve the equations of statistical 
equilibrium — the equations describing how individ- 
ual levels in an atom are populated and depopulated. 
Because the radiation field is no longer Planckian (and 
described by the local electron temperature) the atomic 
populations that satisfy these equations will generally 
differ from their LTE values. In order to solve for the 
populations many atomic processes need to be treated. 
These include: 

1. Photoionization and radiative recombination (i.e., 
bound-free processes) 

2. Low and high temperature dielectronic recombina- 
tion (LTDR and HTDR). In high temperature (or 
classical) dielectronic recombination, recombination 
occurs through d oubly excited autoionizing states 
with large n (e.g.. lBurgesslll964l ). For example, the 
2p nl states in C in that converge on the C IV 2p 
state, and which "recombine" through the decay of 



the 2p electron. In contrast, LTDR usually refers to 
recombination through doubly excited states that lie 
close to, b ut above t he ionization limit (Nussbaumer 
and Storev ll983h . As a consequence, LTDR rates are 
sensitive to the atomic structure — a single energy 
level close to the ionization edge can dominate the 
recombination rate at low temperatures 

3. Bound-bound transitions 

4. Collisional excitation and de-excitation. In hot stars 
this occurs primarily by electrons, but in cooler stars 
collisions with other species are also important. For- 
tunately, the conditions in stellar atmospheres are 
such that the electrons have a Maxwellian velocity 
distribution, and thus we require collision strengths 
averaged over a Maxwellian velocity distribution. 

5. Collisional ionization and collisional recombination 
(i.e., three-body recombination) 

6. Auger ionization by X-rays (and gamma-rays in 
SNe) 

Charge exchange reactions 
Two-photon emission 

In SNe, collisional ionization and excitation with 
non-thermal (high-energy) electrons that are cre- 
ated via gam ma-ray photons and Compt on scatter- 
ing (see, e.g.. lKozma and Fransson 1992 ). 
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When non-LTE is applicable, there is a tremendous 
non-linear coupling between the radiation field and the 
level populations. The radiation field determines the 
electron temperature and level populations which, in 
turn, determine the radiation field. As a consequence 
iterative techniques are used in order to obtain full con- 
sistency between th e radiation fie ld and level popula- 
tions (see review by iHubenvl 120091 ) . 

In Type II SNe, even greater complexities are in- 
troduced. To accurately model Type II SN spectra it 
is important to include advection terms, which couple 
the populations at one time step to those at an earlier 
time step, into the equations of statistical equilibrium. 
The inclusion of advection terms into the rate equations 
helps explain the strength of Ha emissi on in spectra 
of SN 1987A dUtrobin and Chugail 120051; Dessart an d 
Hillierl2010h and SN1999em ( jPessart and Hillierl [2008h . 



4 Atomic Details Matter 

In the construction of model atmospheres the precise 
details of the atomic opacity are generally unimpor- 
tant. However in spectroscopic analyses this is not the 
case. In hot stars we often have only a few lines that 
arise from a given species/ionization stage. In order to 
deduce abundances, it is crucial that our model fully de- 
scribe the formation process of each individual line. In 
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Fig. 1 Illustration of the complex opacity at one location (with T ~ 40, 000 K, iV e ~ 10 15 cm 3 ) for a typical model 
atmosphere. 



LTE this is relatively easy. In non-LTE this can be very 
difficult, since wc need to fully understand the processes 
affecting the populations of the levels involved in the 
transition. This means having all the atomic data rel- 
evant to those levels and, since the population of those 
levels (except in a few special circumstances) depends 
on the population of other levels in the atom, we re- 
quire atomic data for the whole atom. Unfortunately 
the situation can be even more complex than this — in 
some cases line strengths of one species are affected by 
complex interactions with another species. To illustrate 
the complexities, we discuss four examples. 

In WC stars, evolved massive stars deficient in H, 
and showing He, C, and O emission lines in their spec- 
tra, C in A5696 is used as one of the key classification 
diagnostics. In WC4 stars the line is very weak, if not 
absent, while its strength increases (but with scatter) 
as we move along the spectral sequenc e from WC4 to 
WC8 (jTorres. Conti. and Massevlll986f ). Interestingly, 
other C in lines show much less variation in strength 
along the spectral sequence (e.g.. Conti. Massey, and 
Vreux[l99Cj|)~ thus there is something special about the 
formation of C in A5696. 

A simplified Grotrian diagram for the singlet-terms 
of C in is shown in Fig. [2] C in A5696 emission is pro- 
duced by the decay of the 2s 3d *D state to the 2s 3p 1 P° 
state (A — 4.3 x 10 s _1 ). However, photons prefer to 
decay from the 2s3d 1 D to the 2s2p 1 P° state, with 
A = 6.3 x 10 9 s — 1 . In early WC stars, most of the 
decays occur via this route — it is only when this tran- 
sition become s optically th ick, that C in A5696 is driven 
into emission (|Hillierlll989h . To complicate matters fur- 
ther, the strength of C in A5696 is also influenced by 



low-temperature dielectronic recombination. Interest- 
ingly the 2s 3p 1 P° s tate preferentially decays to 2p 2 1 D 
rathe r than 2s 3s 1 S ( Nussbaumerlll971 ; Cardona-Nuned 
Il978l ). Thus A5696 can be in emission, while A8500 
(produced by the decay of 2s Sp 1 ? , the lower level of 
A5696) can remain in absorption — a phenomena seen 
in Of stars (|Ebbets and Wolfjll98lh . 

Another line formation mechanism is continuum flu- 
orescence, which is best illustrated by an example. In 
Civ there is a strong transition at 31 2 A which con- 
nects the ground state, 2s 2 S, to the 3p 2 P° state. As 
we are dealing with a ground state transition, it is usu- 
ally optically thick, and hence photons typically scatter 
many times in this transition before escaping, or before 
being destroyed. However the 3p 2 P° level can also de- 
cay via a transition at AA5801,5812 to the 3s 2 S state. 
The probability of this occurring is low (~ 1/170 per 
scattering) but it does provide a me ans of conve rting 
far UV photons into optical emission ( Hillier 1988| ) . To 
get the strength of the observed optical emission cor- 
rect we have to have a good understanding of the C iv 
model atom as well as accurately model the spectral 
region around 3 12 A — a region which cannot be di- 
rectly observed in O and W-R stars, and which suffers 
heavy line blanketing by iron group elements. Contin- 
uum fluorescence is important in many astrophysical 
objects including W-R stars belonging to the nitrogen 
sequence (WN stars), LBVs, and quasars. 

Perhaps the most famous example of line overlap in- 
fluencing line strengths is seen in nebula spectra. In 
nebula, some O in lines are seen to be unusually strong 
— much stronger than would be predicted by recom- 
bination theory. Moreover, line strengths in individ- 
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ual multiplets do not correspond to those observed in 
the laboratory. The explanation lies in the chance 
coincidence of an Oni line (A303.80) with Hell Lya 
(A303.78). As the He n Lya transition is optically thick, 
photons emitted in the transition are trapped and scat- 
ter many times. During this scattering process there is 
a chance that some of the Lya photons will be absorbed 
by O in. The upper O ill levels have several alternative 
decay routes, some of which lead to the enhanced emis- 
sion in some lines seen at optical and UV wavelengths. 
The process is term ed Bowen Re sonance Fluorescence, 
after its discoverer (|Bowenlll934f ). A more recent and 
detailed discussion of the mechanism is provided bv Qs- 
terbrockll989|). 

In O stars a similar overlap occurs. In this case 
there is a chance overlap of Feiv lines with the Hei 
resonance transition at 584. 33A. The FelV lines remove 
photons from the transition, lowering the population of 
the Is 2p 1 P° level, which in turn affects the strength of 



He I s inglet levels in the optical spectra ([Naiarro et al. 
2006). This in turn affects effective temperature deter- 



minations, since the ratio of He I to He II line strengths 
is used as a temperature diagnostic. Prior to the dis- 
covery of this effect it was known that there were incon- 
sistencies between singlet and triplet He I lines in some 
O stars. The effect was erroneously ascribed (at least in 
this author's opinion) to a problem with the triplet line 
strengths. However, in hindsight it is much more likely 
that the problem lies with the singlet line strengths be- 
cause the lower level of most of the optical diagnostic 
lines is coupled to the Is 2p 1 P° state whose population 
is affected by radiation-transfer effects in the resonance 
transition. 
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5 Atomic Data Requirements 

Astrophysicists are interested in understanding the 
physical processes and properties of astrophysical ob- 
jects (stars, SNe etc). Ideally we would have all the 
required atomic data, and any discrepancies between 
models and observations would only be related to model 
assumptions and the neglect of crucial physical pro- 
cesses. Unfortunately this is not the case — in the real 
world only limited atomic data is available, and it is of 
mixed quality. While great strides have been made im- 
uality and quantity of atomic data, (e.g., 
1987I ; bummer etaH Il993l ; iKurucj HoO^I) , the 



proving 


the 


Seaton 


1987 



Fig. 2 Simplified Grotrian diagram (not to scale) for C III. 



availability of atomic data and its quality must always 
be considered when performing spectroscopic modeling 
to obtain results of astrophysical importance. Given 
limited resources, what are the most crucial data sets 
that are still required? 

In general, the most important elements for hot stars 
are H, He, C, N, O (generally referred to as CNO el- 
ements), and Fe, with Ne, Si, S, and Ar of somewhat 
lesser importance. In SNe the situation is somewhat 
different, with other iron group elements (particularly 
Ni and Co) also being of crucial importance. Even for 
this small subset of species, important atomic data is 
missing. 

One of the most important requirements, but per- 
haps the least appreciated, are accurate energy levels 
and wavelengths. In general, energy levels and wave- 
lengths of sufficient accuracy can generally only be ob- 
tained from observation. As once noted to me, gf values 
accurate to 10% are great but a wavelength accurate to 
10% is (almost) useless. For CNO elements atomic en- 
ergy levels and wavelengths are generally pretty good, 
however additional data are still needed. This is espe- 
cially true of the infrared spectral region. Additionally, 
accurate energy levels are needed for doubly excited 
states that lie near, or above, the ionization limit (e.g. 
C in) since the precise location of these states is impor- 
tant for determining dielectronic recombination rates, 
particularly at low temperatures. Measuring the width 
of lines from these states can also give an estimate of au- 
toionization probabilities, providing an important check 
on theoretical calculations. 

As regards gf values, there has been a proliferation 
of theoretical data which has greatly facilitated the 
advances in spectroscopic analyses. In many cases, 
the calculations assume LS coupling which is gener- 
ally adequate for computing mean opacities and model 
atmosphere structures, but may not be so useful for 



1 Atomic data from Robert Kurucz is available at 
|http://k urucz.harvard.edu/ 
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performing non-LTE abundance studies. Theoretical 
calculations can often provide reliable gf values, but 
for some cases the gf values are not sufficiently accu- 
rate for analyses. As an example we consider the three 
most important decays from the 2p 3p 3 D state in C III. 
These are 

2s2p 3 P° -2p3p 3 D (0.00926, 2.7 x 10 8 s - \ 369.4A), 
2p3s 3 P° -2p3p 3 D (0.166, 1.5 x 10 7 s-\ 6740.6A), 
and 

2s3p 3 P° -2p3p 3 D (0.00286, 4.6 x 10 6 s~\ 1577.4A) 
where the numbers in brackets are the oscillator 
strength, A value, and wavelength, respectively . Usin g 
gf valu es com puted by Nussbaumer and Storevl (|l984T ). 
Hillier f|l987h found that the 1577A line in the WC star 
HD 165763 was much stronger than would be predicted 
on the basis of the observed strength of the 6 740. 6 A 
line. Improved calculations by Peter Storey (private 
communication) indeed showed that the gf value for the 
1577A line was too large, and his revised value (given 
above) gave much better agreement with observation. 
Unfortunately, his calculations also revealed that the 
actual value was quite sensitive to assumptions used in 
the calculations. 

Another example is the gf values of the Fe iv tran- 
sitions which overlap the He 1 resonance transition (see 
Q . The two transitions of most interest are 

n_ „, o j 5 A\ 
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For other elements photoionization cross-sections are 
unavailable. Of particular importance are elements like 
Co and Ni which have high abundances in SNe. At 
present in CMFGEN relatively crud e approximations are 
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3d 4 ( 3 G)4p 2 H° /2 and 
-Jd 4 ( 3 G)4p 4 H° /2 . 



Bell and Kurucz |l995h give gf val ues of 0.00349 and 
0.0288 for the two transitions, while iBecker and Butler 
|l995f ) give 0.00251 and 0.00264. In these cases mea- 
surements are clearly desirable to constrain the strength 
of these lines. 

Photoionization cross sections are now available for 
atoms with even atomic numbers, up to and including 
Fe, p rimarily through the OPACITY an d Iron Projects 
(e.g., Seaton 1987 ; Hummer et al\\199$ ). These cross- 
sections appear to show reasonable agreement with ex- 
periment when available. One possible deficiency in the 
cross-sections is the location of the resonances, which 
are theoretical rather than experimental. Erroneous lo- 
cations can potentially affect recombination rates (par- 
ticularly at low temperatures) and there is also the 
potential problem of interactions between resonances 
and spectral features. The importance of the latter ef- 
fect is unknown. To expedite calculations in CMFGEN 
we generally smooth the photoionization cross-sections 
(typically with a Gaussian of full-width 3000 km s -1 ). 
Originally the smoothing was hard- wired into the cross- 
sections, but for most cross-sections the smoothing is 
now a control parameter. This allows us to easily test 
the influence of the smoothing. As computers become 
faster, and especially for ID models, smoothing is no 
longer a necessity, but rather a computational tool. 



used [ data for Ni 11 is available — iNahar and Bautista 
( 20011 )]. Lines belonging to Sen can also be readily 
identified in Type II SNe spectra, and photoionization 
data for Sc ions is also unavailable. The case of Sc is 
interesting — Sc has relatively low abundance and thus 
it has very little influence on SNe spectra. However it 
does produce readily identifiable features, and match- 
ing such features does provide a quantitative test of 
the spectroscopic model. It also pleases observers (and 
some theorists) , who are quick to point out the absence 
of Sc features from models. 

The lack of collisional data is probably the area of 
most concern. For low lying levels, collisional data is 
generally available, but for higher levels such data is 
usually lacking. In such cases approximate formulae, 
of unknown accuracy, are often used. The formulae 
often depend on the gf value of the transition connect- 
ing the tw o levels, but a s has been pointed out in the 
past (e.g., Mihalasl 19781 ). collisional rates for LS semi- 
forbidden or forbidden transitions can be as large as 
those for non-forbidden transitions. For the vast major- 
ity of applications we need cross-sections averaged over 
a Maxwellian velocity distribution, and it is best tabu- 
lated as a temperature dependent collisional strength. 
D espite recent imp rovemen ts (e.g.. Przvbilla and But- 
ler 20041 ) there are still uncertainties with the hydrogen 
collision strengths. 

Charge exchange reactions are of particular impor- 
tance in non-LTE modeling. In many cases a charge 
exchange reaction with H, rather than the recombina- 
tion rate, determines the ionization balance. An ex- 
cellent example is + + H f± O + H + , which has 
a rate coefficient of order 10 -9 cm 3 s -1 . Many cross- 
sections for charge reaction rates with H and He have 
been computed; a convenient tabulation is that of King- 
don and Ferland (|l996[ ). One problem with tabulated 
cross-sections is that the charge exchange channels are 
not always provided. Such information is necessary as 
we go to high density, since we must include the reverse 
reactions in order to recover LTE. 

Another concern is the paucity of data for charge 
exchange reactions with species other than H and He. 
This paucity is understandable — H and He have the 
greatest effect since they are easily the most abun- 
dant elements in most astrophysical contexts. How- 
ever, there are objects that are H and/or He deficient. 
In WC stars for example, H is absent, and He, C, and 
O have comparable abundances. In SNe a wide range 
of species are present, with C, O, Si, Ne, Fe, Ni, Co 
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exhibiting large mass fractions in some regions. De- 
spite advances in spectr al modeling of SNe (e.g.. Baron, 



gent 20061 ; iKasen and Wooslevll2009l ; iDessart and Hillier 



2010f ). there are still significant uncertainties, and it is 



not currently possible to identify a discrepancy between 
model and observation that might be due to a charge 
exchange reaction. 
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